Fundamental Properties of Eigenstate Wave Functions 
of 'PT— Symmetric Anharmonic Oscillators 
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The PT-symmetric cubic oscillator with Hamiltonian H3 = — ^ d\ + \x 2 +i G x 3 is a paradigmatic 
example of a pseudo-Hermitian, or 'PT-symmetric Hamiltonian with a purely real spectrum when 
endowed with L 2 (R) boundary conditions. Eigcnfunctions of the stationary Schrodinger equation 
H3 'tpn(x) = E n ip n (x) are manifestly complex, while the energy eigenvalues E„ are real and positive. 
Although Hi, does not commute with the parity operator, we find that for a natural choice of the 
global complex phase of the eigenstate wave function, the real and imaginary parts of the eigenfunc- 
tions [i.e., Keip n (x) and lmtp n (x)] are eigenstates of parity. Both the real and imaginary parts of 
the eigenfunctions are found to have an infinite number of zeros, even for the ground state and even 
for infinitesimal G, but the real and imaginary parts of ip„(x) never vanish simultaneously when 
G > 0. Furthermore, we find that the eigenfunctions are "concentrated" in an "allowed" region 
where the energy E n is larger than the complex modulus of the complex potential. PT-symmetric 
Hamiltonians constitute natural generalizations of Hermitian Hamiltonians as time-evolution op- 
erators. Our results suggest that PT-symmetric (pseudo-Hermitian) time evolution can naturally 
be interpreted as time evolution in a situation where manifestly complex "gain" and "loss" terms 
mutually compensate and lead to the real energy eigenvalues. 
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I. INTRODUCTION 

Quantum mechanics promotes a number of physical 
quantities, which are purely real in classical physics, into 
the complex domain. An example is the momentum 
which is a real, measurable physical quantity in classical 
physics, but its quantum mechanical equivalent, the mo- 
mentum operator p = — i V, contains the imaginary unit. 
A quantum mechanical wave function ip(r) is manifestly 
complex. Plane-wave eigenstates of the free Schrodinger 
Hamiltonian Hq = — V 2 /(2m) also are eigenstates of 
the momentum operator p = — iV and are of the form 
exp(i k ■ r). A purely real eigenfunction such as cos(fc • r) 
does not describe a particle with definite momentum. 

Yet, under special circumstances, it is possible to con- 
struct eigenfunctions of operators which are purely real. 
Indeed, it is known [1, 2] that the eigenfunctions of the 
harmonic oscillator can be chosen as purely real func- 
tions, and furthermore, eigenfunctions of the quartic an- 
harmonic one-dimensional oscillator 



H4 — — — d 2 
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gx 



9>0, 



(1) 



can also be chosen as purely real. Here, we assume that 
H4 is endowed with £ 2 (1R) boundary conditions. This 
implies that the wave function vanishes at infinity, i.e., 
— > for x — > ±00. 
It has recently been observed [3-7] that the Hamilto- 
nian 



Ho 



1 



d 2 . 



1 



iGx 3 , G>0. 



(2) 



has a real spectrum when endowed with L 2 (R) bound- 
ary conditions. The Hamiltonian H3 commutes with the 



combined parity-time-reversal operator, \P T, H3] = 0. 
Here, V is the parity x — > — x, and T is the time rever- 
sal, which for H3 amounts to the replacement i — > — i. 
Operators with these properties have been analyzed re- 
cently [3, 8] and have been advocated as viable alterna- 
tives to Hermitian operators. The TT-symmetry prop- 
erty implies that 



H 3 = VTH 3 rV = VH+V. 



(3) 



where H 3 is the Hermitian adjoint. Here, we have used 
the relation V 2 = T 2 = 1 and thus (TT) _1 = TV. 

If a Hamiltonian H fulfills a relation of the type H = 
j/ -1 H + 77, then H is said to be pseudo-Hermitian [9]. 
Here, following Pauli [9], we only require the metric op- 
erator 77 to be Hermitian. The Hermitian adjoint H + is 
calculated as H + = THT, and thus, TT-symmetry is 
identified as a special case of pscudo-Hcrmiticity (rj = V). 
One can generically show [9] that the scalar product 



dx tjj* (t, x)V(j)(t, x) 



(4) 



is conserved under time evolution if both ip and <fi ful- 
fill the time-dependent Schrodinger equation idtip{t) = 
H3tjj(t), and idt<fi(t) = H^(j>{t). However, the in- 
tegrand in Eq. (4) is manifestly "nonlocal" because 
ip*(t, x)V(j){t, x) = tp*(t,x)(j)(t 1 —x); it depends on func- 
tions evaluate at x and —x. 

In the prolific literature on PT-symmetric quantum 
mechanics (see Refs. [3, 5, 6, 8-18] for an incomplete list), 
the spectrum of PT-symmetric Hamiltonians has been 
analyzed in detail. However, much less attention has 
been devoted to the wave functions corresponding to the 
eigenstates. An analysis of the wave functions appears in- 
dicated because a number of interesting field-theoretical 
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FIG. 1: (Color online.) The ground-state energy of the anharmonic oscillator Hamiltonian with imaginary cubic 
perturbation, as given in Eq. (2), is plotted in the perturbative (weak-coupling) regime in Fig. (a). The solid curve 
represents the results from eighth-order perturbation theory (terms up to the order G 8 ). In Fig. (b), the same 
eigenvalue is plotted in the strong-coupling regime G > 0. The small dots represent numerical values. 



model theories and a streamlined description of phe- 
nomenologically important so-called PT-symmetric wave 
guides rely on PT-symmetric quantum mechanics and 
field theory (a brief overview is provided in Appendix A) . 
Moreover, the concept of a PT-symmetric Hamiltonian 
has recently been instrumental the Bender-Wu formu- 
las [19-21] for odd anharmonic oscillators [5, 6]. It is 
self-evident that an improved intuitive understanding of 
the physics involved in PT-symmetric models requires 
further understanding of the wave functions. 

With reference to the wave functions, the following 
question is of interest: The Hamiltonian H 3 given in 
Eq. (2) involves a manifestly complex potential, which 
we denote as W(x), 



1 



W{x) = -x 2 +iGx 3 = F(x)e i 



arg(V(») 



V{x) = \W(x)\ = yf\ X* + G 2 X 6 



(5a) 



(5b) 



The complex modulus V(x) = \ W(x)\ tends to infinity as 
x — > ±oo. For purely real potentials like the "confining" 
quartic oscillator given in Eq. (1), intuition suggests that 
the "bulk" of the probability density of the eigenstatc 
wave function should be concentrated in the "classically 
allowed" region, i.e., in the region where the eigenenergy 
E is greater than the potential, i.e., E > V(x) [where 
V(x) £ R]. The generalization of this concept to PT- 
symmetric quantum mechanics needs to be clarified. For 
a manifestly complex potential, the condition E > W(x) 
[with W(x) s <D] cannot be applied because the complex 
numbers are not ordered. 

The quartic anharmonic oscillator Hamiltonian given 
in Eq. (1) commutes with parity. By contrast, the eigen- 
state wave functions of a PT-symmetric Hamiltonian do 
not need to be eigenstates of parity, because the par- 
ity operator does not commute with the PT-symmetric 



Hamiltonian. These observations raise a number of ob- 
vious pertinent questions which we aim to address. We 
proceed by recalling a few basic facts about eigenvalue 
perturbation theory in Sec. II, analyzing the parity of 
eigenstates in Sec. Ill, and continue with a visualization 
of the manifestly complex PT-symmetric eigenstates in 
Sec. IV. This leads to a generalization of the "classi- 
cally allowed" region to the PT-symmetric domain. We 
conclude the discussion in Sec. V. 



II. ASYMPTOTICS OF IMAGINARY CUBIC 
PERTURBATION 

In order to fix ideas, we first recall a few basic facts 
about eigenvalue perturbation theory and the imaginary 
cubic perturbation. One can split the Hamiltonian given 
in Eq. (2) into an unperturbed harmonic-oscillator term 
iT°) = -\dl + \ x 2 and a perturbation = iGx 3 . 
Starting from the unperturbed harmonic oscillator eigen- 



values H^\n) = en \n), where e„ = n + 5, one can 
develop perturbation theory for the eigenvalues e„ of 
H 3 = + either using the classical Rayleigh- 

Schrodinger approach [2] or using a complex contour in- 
tegration of the logarithm of the wave function (Riccati 
equation, sec Rcfs. [22, 23]). The first nonvanishing term 
in perturbation theory is of second order and reads as 
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FIG. 2: (Color online.) Figure (a) displays the ground-state wave function of the imaginary cubic Hamiltonian (2) 
for G = 1.0. This wave function is manifestly complex. The real part Reipo{%) (even under parity) is plotted using 
solid lines, and the dashed curve corresponds to the parity-odd imaginary part Im^o(i). For the first excited state 
(still, G — 1.0), the real part is odd, while the imaginary part is even under parity [see Fig. (b)]. The second 
excited state [Fig. (c)] has an even real part, while its imaginary part is odd. The global complex phase of the wave 
function is chosen so that the real part Re ip n =o, 1,2(2;) of the wave functions has the same qualitative behavior as 
the eigenstate wave function of the quartic oscillator displayed in Fig. 4(a). 



so that e„ ps e„ + e„ . There is an observation to be 
made: Unlike for any second-order perturbation involv- 
ing a Hermitian operator, the second-order term here is 
positive and shifts the ground-state energy level upward 
instead of downward. Through fourth order, the result 
reads as 

1 G 2 

e„-n+- + — (30n 2 + 30n+ll) (7) 

2 8 

~ S ° 4 ( 94n3 + 141n2 + 109n + 21 ) + £>(G 6 )) , 

where n is the principal quantum number. For the 
ground state (n = 0), this expression evaluates to 
1/2 + llG 2 /8 - 465G 4 /32, while going to eighth or- 
der we obtain 1/2 + 11 G 2 /8 - 465 G 4 + 39709 G 6 /128 - 
19250805 G 8 /2048 which is plotted against numerical val- 
ues of the ground-state energy in Fig. 1(a). The posi- 
tive curvature of the ground-state energy in the weak- 
coupling regime is clearly visible. 

We can also look at the strong-coupling asymptotic as 
G — > 00, for which wc employ the "poor man's scaling" 
x — > G~ 2 / b x. The Hamiltonian H% is thus replaced by 
i?3, which has the same spectrum as H3, 

H' z = Gr.{- l -dl + ^ + \G~'i X ^ . (8) 

For large G, the eigenvalues e n of therefore go as 

~ d e n , G -> 00 . (9) 

Here, the e n are the energy eigenvalues of the Hamilto- 
nian H3, with 

# 3 = -^ 2 + ia; 3 ^2"l(-G) 2 + ix 3 ) . (10) 

Indeed, the 'PT-symmetric Hamiltonian H3 has its own 
set of eigenvalues, which we label as e n . In the last re- 
placement step in Eq. (10), we have done the scaling 



transformation x — > x. This transformation allows 
us to connect the strong-coupling asymptotics with the 
literature, notably, with Refs. [3] and [4, 7]. In particu- 
lar, we have e n = 2r'~> E n where the E n are the energies 
of the Hamiltonian — <9 2 + ix 3 . 

III. PARITY OF EIGENSTATES 

We have determined, to high numerical accuracy, the 
manifestly complex wave functions of the ground state 
and the first three excited eigenstates of the imaginary 
cubic oscillator (see Fig. 2). Because the parity opera- 
tor V does not commute with the Hamiltonian H3, the 
eigenstates of H 3 are not eigenstates of parity. Further- 
more, because the potential is manifestly complex, so 
are the wave functions. Yet, numerical evidence drawn 
from Fig. 2 suggests that individually, both real as well as 
imaginary part of the eigenstate wave functions are eigen- 
states of parity. Eigenfunctions of a quantum Hamilto- 
nian are only defined up to a global phase factor exp) (iip) 
(and a dilation factor, in the sense of "rays" in Hilbcrt 
space)). The above statement therefore needs to be qual- 
ified: It holds if we perform the numerical diagonalization 
of the Hamiltonian (2) in a basis of harmonic-oscillator 
eigenfunctions. 

We can formally split the Hamiltonian H3 into a "real 
part" and a "imaginary part" as follows, 

RcH 3 = - \d 2 x + ix 2 , (11a) 

Imff 3 = iGa: 3 . (lib) 

If we also split the eigenstate wave function ipn(x) into 
real and imaginary parts, 

ip n (x) = Re i/) n (x) +ilm ip n (x) , (12) 

then it is rather easy to show, based on the decompo- 
sition (11), that if Hetp n (x) is even under parity and 
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ip n {x) is an eigenstate of if 3 with real eigenvalue of e„, 
then Im^ n (x) has to be parity-odd, and vice versa. Also, 
if Retpn(x) is odd under parity and "ip n {x) and the corre- 
sponding eigenvalue e„, is real, then Imi/) n (i) has to be 
even under parity. 

In Ref. [24] , it has been observed that since VT com- 
mutes with the Hamiltonian, the eigenfunctions of H3 
also have to be eigenfunctions of the VT operator. The 
precise eigenvalue of VT may, however, depend on the 
phase assigned to ip n (x) because T is an antilinear oper- 
ator. Let us first investigate the phase conventions used 
in Fig. 2, where the real and imaginary parts of the wave 
function are, alternatingly, even and odd under parity 
as we proceed to higher excited states. The appropriate 
eigenvalues are thus 

VT^ n {x)^iC n {x) = {-lT^{x). (13) 

In this convention, the real and imaginary parts of the 
eigenfunctions are alternatingly even and odd, as we 
would expect for anharmonic oscillators [1, 2]. 

However, we may simply adjust the phase of the eigen- 
statcs as follows, 

VT[(-i) n M*)] = i" C(z) - [(-i)" 1>(x)] , (14) 

and so the wave functions [(— i) n ip n (x)], which we might 
just as well denote by ip n (x), are eigenstates of H3 
and eigenstates of 'PTwith eigenvalue +1 [cf. Eq. (5) of 
Ref. [24]]. With reference to Fig. 5, let us now define the 
functions $ n (%) = ( — i)™^™^) as the new eigenfunctions. 
In the space of the *]/„, the conserved scalar product (4) 
then becomes 

<*„|*m>. = / dx**(x)VV m (x) 

= jdx (VT*n) (*) * m (z) 

= J dx^„ {x)^ m (x) , (15) 

where the last line is without complex conjugation. This 
latter formula shows that within the space of eigenstates, 
the conserved scalar product for the 'PT-symmctric imag- 
inary cubic perturbation is equal to the so-called gener- 
alized inner product for complex-scaled Hamiltonians as 
defined in Eq. (2.4.2) of Ref. [25], which avoids the com- 
plex conjugation of the first argument. Also, the inte- 
grand in Eq. (15) is "local" in the sense that it depends 
only on function values of argument x, not —x. 

We thus identify the 'PT-symmetric scalar product as 
the generalized inner product for complex-scaled Hamil- 
tonians (which otherwise give rise to resonances (see 
Refs. [5, 6, 25]). For a "real" cubic perturbation G x 3 in- 
stead of an imaginary cubic perturbation \G x 3 , the gen- 
eralized inner product describes the orthogonality prop- 
erties of resonance wave functions. [4, 25], corresponding 
to the manifestly complex resonance energies [5, 6]. 
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FIG. 3: (Color online.) Modulus-phase plot of the 
harmonic+cubic potential given in Eq. (5a). The modulus 

V(x) — \W(x)\ — ^ j x 4 + G 2 x 6 leads to a confining mech- 
anism for x — » ±00. The value of G in the plot is G = 1.0. 
The shaded area displays the complex phase of the potential 
and covers the interval (— n/2, 7r/2). 



IV. VISUALIZATION OF 'PT— SYMMETRIC 
EIGENSTATES 

Because the energy eigenvalues of the imaginary cubic 
oscillator (2) are real, we can associate the imaginary cu- 
bic perturbation with the real, "confining" quartic poten- 
tial of the quartic oscillator given in Eq. (1). In the latter 
case, the wave functions are known to be "confined" to 
the "classically allowed" region where E > V(x). Fur- 
thermore, from the investigations [5, 6, 10], we know 
that a dispersion relation connects the energy levels of 
the imaginary cubic perturbation iG x 3 to the "real" cu- 
bic perturbation Gx 3 . However, the eigenenergies of the 
imaginary cubic potential are real, whereas the resonance 
and antiresonance energies of the real cubic perturbation 
are complex [4]. At some risk to oversimplification, we 
can say that the role of the stable quartic perturbation 
+g x 4 in Eq. (1) is taken over by the imaginary cubic per- 
turbation i G x 3 , whereas the role of the unstable quartic 
perturbation —g x 4 is played by the real cubic perturba- 
tion G x 3 (here, we assume g > and G > 0). There is no 
direct and obvious visualization available for the "clas- 
sically allowed" region in an imaginary cubic potential. 
One may therefore ask in which sense the imaginary cu- 
bic perturbation "confines" the eigenstate wave functions 
to a "classically allowed" region of space, in particular, 
because the set of complex numbers (D cannot be ordered. 

An intuitive understanding can be obtained if we in- 
terpret the potential in terms of a complex modulus and 
a phase, which according to Eq. (5a) reads as 

W(x) = i x 2 + i G x 3 = V(x) c 1 ai s(v(x)) ^ ( 16a ) 

-tt/2 < axg(W(x)) < tt/2 . (16b) 

In Fig. 3, we plot the modulus of the complex potential 
and its complex phase. 

The eigenstate wave functions are plotted in Figs. 4, 5 
and 6, with the idea of writing the complex eigenstate 
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(a) (b) 

FIG. 4: (Color.) In Fig. (a), we plot the wave function of the ground and the first two excited states of the quartic 
oscillator. The "stable" quartic perturbation perturbs the potential in Fig. (a) and leads to a "confinement" of the 
wave function to the classically allowed region E > V(x) = |s 2 + g x 4 . The eigenenergies of the imaginary cubic 
perturbation are real and allow us to plot the complex "PT-symmetric eigenfunctions in Fig. (b). For the latter, we use 
a modulus-phase-plot in analogy with the potential W(x) displayed in Fig. 3. The complex phase of the wave function 
is displayed in the shaded region. The complex phase 9 — 6(x) in the decomposition ip n (x) = ^^(a;)! exp[i#(x-)] 
covers the interval (— 7r,7r). The red curves denote the complex modulus square of the wave functions [Figs, (a) 
and (b)], while the blue curves display the complex phase of the wave functions. The values of g and G in the plots 
read as 1.0 and 1.0, but it is useful to note that the qualitative behavior of the curves does not depend on the concrete 
value of g. For the PT-symmetric potential [Fig. (b)], the quantity [^(x)! 2 does not have an interpretation as a 
conserved probability density, but together with the phase information, it specifies the wave function uniquely. The 
qualitative features ( "humps" ) of the quartic oscillator are still present in the complex ( "PT-symmetric" ) domain, 
however the zeros of the wave functions merely become local minima. See also the following Fig. 5. 



wave functions as 

Mx) = \Mx)\e iarsi ^ {x)) ■ (17) 

It is known that for Hcrmitian operators, the ground- 
state wave function always has maximum symmetry. Fur- 
thermore, for the "stable" quartic perturbation, the num- 
ber of zeros of the wave function is equal to the princi- 
pal quantum number. This is illustrated in Fig. 4(a). 
For the "unstable" , "real" cubic perturbation, it has al- 
ready been observed in Ref. [4] that the resonance state 
wave function has no zeros when considered as a com- 
plex variable (see Fig. 4 of Ref. [4]). However, a priori, 
this statement cannot automatically be generalized to the 
"stable", imaginary cubic perturbation. In Fig. 4(b), the 
plot of the wave function square If/ViXaOl 2 suggests that 
indeed, the same statement holds for the eigenstate wave 
functions of the imaginary cubic perturbation; they have 
no complex zeros. We have used a numerical value of 
G = 1.0 in Fig. 4; the ordinate axis covers the range 
< V(x) < Vo — 7. In Fig. 6, we illustrate that the "con- 
finement" mechanism holds for all values in the range 
< G < 1. The wave function squares are plotted in ar- 



bitrary units; in practical calculations, one normalizes to 
/ dxlV^s;)! 2 = I- The complex phase of the wave func- 
tion, which covers the shaded areas in Fig. 4(b), revolves 
in the complex plane and transits a number of Riemann 
sheets, i.e., it jumps from —it to +7r several times in the 
interval —3 < x < 3. 

A Wentzel-Kramers-Brioulin (WKB) analysis shows 
that the phase of the wave function has to behave, asymp- 
totically, as a,Tg(tp n (x)) ~ —x 5 / 2 for x — > +oo and as 
arg( - 0„(x)) ~ x 5 / 2 for x — > — 00, and we have verified this 
asymptotic behavior against our numerically determined 
wave functions. Every time the complex phase of the 
eigenstate wave function a,rg(tp n (x)) (modulo tt) attains 
zero, the imaginary part of the wave function vanishes, 
and whenever &Tg(T/j n (x)) mod n = ir/2, the real part of 
the wave function vanishes. So, we conclude that, even 
for an infinitesimally small coupling G in the imaginary 
cubic perturbation i G x 3 , the real and imaginary parts of 
the wave function, individually, have an infinite number 
of zeros. This is somewhat surprising. E.g., the above 
considerations imply, among other things, that the num- 
ber of complex zeros of the first-excited-state wave func- 
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FIG. 5: (Color.) In Fig. (a), we plot the probability density p = |?/>(a;)| 2 of the quartic oscillator's ground state and 
first two excited states, in a potential V(x) = ^x 2 + g x 4 with g = 1.0, but in a different way as compared to Fig. 4. 
Although the wave functions of the quartic potential are purely real, we use a complex modulus-phase plot for the 
real wave functions in Fig. (a). A sign change then corresponds to a jump in the complex phase from zero (for 
a positive real number) to —n (for a positive negative number). In Fig. (b), we give a modulus- phase plot of the 
eigenstate wave functions of the imaginary cubic perturbation, but this time we shift the global phase of the wave 
function as compared to Fig. 4(b) and Fig. 2, by multiplying it with a phase factor (— i) n [see also the discussion 
surrounding Eq. (15)]. Here, the complex phase 9 — 6(x) enters the decomposition (— i) n if) n (x) = \ip n (x)\ exp[i#(x)] 
and covers the interval (— 7T,7r). 



tion of the imaginary cubic perturbation is discontinuous 
at G = 0. Namely, for G = 0, we have one complex zero 
(because the Hamiltonian is equal to the harmonic oscil- 
lator). For novanishing G, the total number of complex 
zeros of the wave function vanishes, while both real and 
imaginary parts of the wave function (for nonvanishing 
G) have an infinite number of zeros. 



V. CONCLUSIONS 

Let us add some concluding remarks. It has been 
stressed in the literature that the scalar product {ip\V\(f))* 
as defined in Eq. (4) is not positive definite. This 
has been used as an argument against the viability of 
■PT-symmctric Hamiltonians for the description of nat- 
ural phenomena. However, one may counter argue that 
the same problem persists with regard to the relativis- 
tic Klein-Gordon equation where the time-like compo- 
nent of the current can become negative (see Chap. 2 of 
Ref. [26]). The Klein-Gordon equation is supposed de- 
scribes a charged scalar field like the charged component 
of the Higgs (doublet) field [27] (the latter is usually as- 
sumed to vanish under a gauge transformation, and the 
remaining neutral component of the Higgs doublet is ex- 
panded about its vacuum expectation value). Strictly 



speaking, one has to reinterpret the zero component of 
the conserved Noether current as a charge density, not a 
probability density. Analogously, in the context of neu- 
trino physics, an interpretation of the (non-positive defi- 
nite) timelike component of conserved Noether current of 
the pseudo-Hcrmitian, generalized Dirac interaction den- 
sity as a non-positive definite "weak-interaction density" 
(of the neutrino) has recently been proposed [28] . 

No physical system is known which is directly de- 
scribed by an imaginary cubic perturbation, and so 
the question of how to interpret the conserved density 
P*(x) = ip(x) + 'Pip(x) = ip{x) + ^i{—x) is perhaps not as 
immediate as in the other cases. In Eq. (15), we show 
that the integrand of the conserved current can be rewrit- 
ten in terms of a "local" expression, containing only func- 
tion evaluations at x, and interpreted in terms of a gen- 
eralized inner product [25], which refrains from any com- 
plex conjugation under the integral sign [see Eq. (15)]. 

Incidentally, the reality of the eigenvalues of the imagi- 
nary cubic perturbation is somewhat tied to the discrete- 
ness of the spectrum: Namely, for G = 0, one starts from 
the unperturbed harmonic oscillator. In general, for a 
■pX-symmetric Hamiltonian, one can only show that if 
e„ is an eigenvalue, so is the complex conjugate e* . How- 
ever, if we assume that all eigenvalues are smooth func- 
tions of G in the limit G — > 0, without level crossings or 
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FIG. 6: (Color.) Illustration of the confinement mechanism 
for the imaginary cubic potential described by the Hamilto- 
nian (2), for the ground and the first excited state. The bulk 
of the modulus square of the wave function is centered in the 
"allowed" region where the (real rather than complex) energy 
E > V(x) = V(G,x) = \W(G,x)\. The potential is plotted 
in green, the complex moduli of the wave functions are plot- 
ted in green. The ground state wave function has a modulus 
square \i/jo(x)\ 2 = \4>o(G, x) 2 as a function of G and x. As G 
increases, the bound-state energy (which is equal to the base 
line of the wave function curve at any given value of G) in- 
creases, and the modulus of the potential forms a more narrow 
trough to which the ground-state wave function is confined. 
The same is true for the first excited state. The central mini- 
mum of the complex modulus square of the first-excited state 
wave function is clearly visible. 



splittings or other exotic phenomena, then the reality of 
the eigenvalues for nonvanishing G follows immediately. 
The latter situation is commonly referred to as "broken" 
PT-symmetry [24]. 

For the imaginary cubic perturbation, a "Hcrmitizing" 
transformation, which leads to an equivalent Hermitian 
Hamiltonian has been discussed in [24, 29, 30]. How- 
ever, our discussion in Appendix A suggests that non- 
Hermitian time evolution operators can be useful for 
the description of physics phenomena in specific cases 
where absorption and emission phenomena have to be 
described. An example is given by PT-symmetric wave 
guides [31-33]. Let us interpret terms in the Hamilto- 
nian with a negative imaginary part as "absorptive" or 
"loss" terms and terms with a positive imaginary part as 
"emissive" or "gain" terms. Terms of this nature are of- 
ten used in laser physics (interactions of few-level atomic 
systems with light fields [34]). Our discussion in Sec. IV 
suggests that PT-symmetric Hamiltonians can naturally 
be interpreted as time-evolution operators which describe 



physical systems where the "emission and absorption pro- 
cesses are in equilibrium" (figuratively speaking) or "mu- 
tually compensate each other" , and the eigenenergies of 
the system therefore are purely real. A typical behav- 
ior of such eigenenergies and corresponding wave func- 
tions is illustrated in Figs. 2 — 6. We find that the real 
and imaginary parts of the complex eigenstate wave func- 
tions are still eigenstates of parity, individually, even if 
the 'PT-symmetric Hamiltonian (2) does not commute 
with parity. The "confining" nature of the complex mod- 
ulus of the cubic anharmonic oscillator is illustrated in 
Fig. 3. The wave functions are concentrated in the "al- 
lowed" region with E > V(x) = |W(x)|, as illustrated in 
Figs. 4—6. 

Our considerations, reported in Figs. 4 — 6, suggest 
that the eigenstate wave functions of PT-symmetric 
Hamiltonians are not as "exotic" as one might other- 
wise imagine. An interesting observation, described in 
Sec. IV, is as follows: In a PT-symmetric case, the num- 
ber of zeros of the wave function cannot be used to enu- 
merate the eigenstates. Both the real as well as the imag- 
inary part have an infinite number of zeros, individually, 
but they never vanish simultaneously. The modulus of 
the potential, which tends to infinity as \x\ — > oo, is re- 
sponsible for the confinement of the wave function of the 
imaginary cubic Hamiltonian to a "classically allowed re- 
gion" E > V(x) = |W(a;)| (where W may be complex). 
Finally, we note that the field of PT-symmetric Hamil- 
tonians gives rise to a (still unanswered) number of the- 
oretical questions, even if a number of the ideas used in 
the field have originally been formulated 70 years ago [9] . 
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Appendix A: Background Information 

At some risk to oversimplification, one can point out 
that the scientific world would be a much more interest- 
ing place if the restriction on the Hcrmiticity of time- 
evolution operators could be relaxed somehow. E.g., a 
complex refractive index of a material describes an ab- 
sorption process in which a light wave loses intensity as 
it travels through a material. In 1943, Pauli introduced 
the idea of pseudo-Hermiticity in which the restrictions 
for the Hamiltonians is relaxed [9]. Pauli required that 
a pseudo-Hermitian operator H fulfills a relation of the 
form H = rf 1 H + r\. If r\ = 1, then the Hamiltonian H 
is Hermitian. Along the same vein, in 1998 Bender and 
Boettchcr proposed PT-symmctry [3]. Existing strong 
numerical evidence suggests that the spectra of many 
typical PT-symmetric Hamiltonians are real or come in 
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complex conjugate pairs (the latter being referred to as 
"broken" PT-symmetry). Hcrmiticity does not imply 
PT-symmetry. The PT-symmetric Hamiltonians consti- 
tute a new class of time-evolution operators. 

Recent years have seen both PT-symmetry and 
pscudo-Hermiticity gain a foothold in quantum mechan- 
ics and field theories, including the following four areas, 
(i) Several quantum field-theoretical models, which had 
been deemed problematic under the restrictions imposed 
by Hermiticity, were shown to be consistent under a PT- 
symmetric interpretation. Notably, the ghost state in the 
Lee model was shown to have a positive norm when suit- 
ably reinterpreted [35]. (ii) In cosmological models, it is 
sometimes useful to write down a scalar field with nega- 
tive kinetic energy. Unfortunately, this method leads to 
inherently unstable theories. Andrianov et al. [36] have 
studied cosmological models coupling two fields, one of 
which has a complex potential while both have a positive 
kinetic term. This model is described by a PTsymmetric 
Lagrangian. Furthermore, this model may help resolve 
a number of paradoxes in the evolution of the Universe 
from the big bang to the "big rip" (see Rcf. [36]). (iii) 
The index of refraction is used to describe the propa- 
gation of light through a medium. When the medium 
is opaque, light is absorbed and the index of refraction 
becomes complex. As such, the evolution equations of 
certain waveguides are PT-symmetric Schrodinger equa- 



tions. One such example is a double channel waveguide, 
with one "loss" and one "gain" channel [32] . The "time" 
direction is the direction of propagation. PT-symmetric 
honeycomb lattices [31] have also been studied. By intro- 
ducing an alternating gain-loss structure on a hexagonal 
optical lattice, the dispersion relation can be tailored to 
be equivalent to that of tachyons [31, 37-41]. (iv) Cho- 
dos, Hauser and Kostelecky [42] have proposed a supcr- 
luminal (tachyonic) model of the neutrino which leads to 
a natural suppression of left-handed neutrino, and right- 
handed anti-neutrino states. The PT-symmetric char- 
acter of the model has recently been clarified [43], and 
a number of "sum rules" have been found which lead 
to tremendous simplifications in the calculation of the 
tachyonic propagator [28, 44], and may even offer natu- 
ral explanations for hitherto unexplained paradoxes con- 
nected to neutrino physics [28]. It should be stressed that 
current experimental data neither excludes nor confirms 
neutrino propagation exceeding the speed of light [28, 45] . 
Final clarification can only be obtained once more precise 
measurements of the neutrino propagation velocity and 
the beta decay endpoint spectrum are performed. The 
best estimates obtained from all recent measurements 
of the neutrino mass square using beta decay endpoints 
have been negative (neutrino mass square mf, < 0), while 
still being consistent with zero within experimental error 
(for a list of experimental values, see Ref. [46]). 
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